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A complete explanation of the synthesis of metal-catalyst nanoparticles, and the subsequent nu- 
cleation and growth of bundles of C-SWNTs is introduced using a novel model. It is shown that 
the synthesis process leads to the formation of a liquid layer supersaturated in carbon surrounding 
each metallic-catalyst nanoparticle. The onset of a solutal Benard-Marangoni instability and the 
subsequent formation of patterns of hexagonal convection cells in the liquid layer is predicted and 
quantified by linear and weakly nonlinear analyses. The nucleation and growth of a C-SWNT at 
the center of convection cell is explained. 
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This Letter presents a model which describes the mech- 
anism for the nucleation and growth of carbon single- wall 
nanotubes (C-SWNTs). Several models have been pro- 
posed in the literature Q, 0, 0, 0] ; they have succeeded 
in describing the general scenario leading to growth of 
C-SWNTs. Our model, however, is able to explain spe- 
cific and important aspects such as the nucleation, the 
individual structure of C-SWNTs and their diameter. In 
brief, our model explains how: 1. the synthesis conditions 
for C-SWNTs lead to the formation of carbide-metal 
nanoparticles surrounded by a layer of liquid supersat- 
urated in carbon; 2. a Benard-Marangoni instability is 
generated in this liquid layer; 3. the Benard-Marangoni 
instability leads to the formation of hexagonal convection 
cells, at the center of which nanotubes nucleate; 4. the 
C-SWNTs grow by a root-growth mechanism at least ini- 
tially and 5. the BMI offers a simple kinetic explanation 
for the diameter of the C-SWNTs. 

Recently, a review of the nucleation and growth of 
C-SWNTs suggested that an interface instability is 
responsible for the nucleation of C-SWNTs. Leverag- 
ing this idea, we develop a model based on the Benard- 
Marangoni instability (BMI), called the BMI model. In 
this Letter, we show that nucleation of C-SWNTs indeed 
results from the onset of the BMI. For concreteness, we 
study the common case of a long bundle growing on one 
side of a large (5nm to 20nm) nanoparticle (NP) and only 
consider an iron catalyst, which is known to result in syn- 
thesis of C-SWNTs The treatment remains valid for 
any catalyst. 

We first describe the formation of the nanoparticle and 
the liquid layer surrounding it, by analysing the iron and 
carbon (Fe-C) phase diagram (FigUJ. In gas-phase pro- 
cesses, the synthesis begins by the creation of a mixture of 
Fe-C vapors in an inert gas. The Fe-C mixture is initially 
at a high temperature (5000°C) (point P0 in FigQJ and 
contains a fraction of about 2%-3% of iron atoms with re- 
spect to carbon atoms. The vapor mixture is then cooled 



very rapidly (« 10 6 °C/s) . As the temperature decreases 
to around 2300°C (PO to PI in Fig.l), the vapor mixture 
condenses to form nanometric droplets of an iron-carbon 
solution having a concentration of about 25%-30%at. If 
the cooling rate remains high, the iron-carbon droplets 
undergo a non-equilibrium solidification process down to 
around 1500°C, maintaining droplets in a supersaturated 
state (PI to P2 in FigQ. Subsequently, droplets experi- 
ence a rapid segregation process |3j so as to recover their 
equilibrium state. At this point, the carbon concentra- 
tion in the droplet is higher than the equilibrium liquidus 
value. As a consequence, carbon is expelled radially to- 
ward the surface of the droplet, until the composition 
of the droplet center reaches the liquidus equilibrium. 
From this point, the core of the droplet follows a differ- 
ent path in the phase diagram than the layer surrounding 
this core. 

We now discuss the evolution of the droplet core. This 
core solidifies at the eutectic temperature, T eut (P2 to 
P3 in FigQ, by a eutectic transformation that results in 
the formation of a carbide (Fe^C) phase and of a 7 — Fe 
phase. In agreement with this analysis, the presence of 
a thin carbide layer has already been observed P, by 
transmission electron microscopy at the surface of the 
core. The analysis is further supported by Saito 0] and 
Alvarez et al. Q, who showed that the synthesis tem- 
perature of C-SWNTs and their growth are correlated to 
the catalyst-carbon eutectic temperature. 

In the outer layer of the droplet, the expelled carbon 
will not have enough time to nucleate into small graphitic 
flakes. Carbon remains dissolved in atomic form in the 
liquid metal, resulting in the formation of a nanometric 
liquid layer surrounding the solid core. This liquid layer 
is supersaturated in carbon with a concentration of 
as much as 40%-50%at.(P2 to P4 in FigQ}. We note that 
the liquidphase still exists few hundreds of degrees be- 
low T eut @, 0] since a solution having a non-equilibrium 
carbon concentration (> 25%at. for Fe-C, Ni-C and Co- 
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FIG. 1: The overall scenario for growth of C-SWNTs repre- 
sented in the Fe-C phase diagram 

C binary systems) has a lower melting temperature than 
the eutectic temperature of the system at equilibrium. 
We can thus define a non-equilibrium eutectic tempera- 
ture for the supersaturated solution. The melting point 
of this solution can vary from 800°C to 1000°C and will 
depend on the carbon concentration, a higher concen- 
tration leading to a lower melting point. Therefore, the 
surface liquid layer can exist from P4 to P5 in the phase 
diagram FigQ] This is also supported by Gorbunov et 
al. [2| who proposed that the catalytic nanoparticles in- 
volved in the C-SWNTs synthesis are molten. 

We now focus on a single NP, and explain a mecha- 
nism for the nucleation and growth of nanotubes, driven 
by a convective instability. A NP is composed of a liquid 
layer supersaturated in carbon surrounding a solid core 
composed of a carbide phase and of a 7 — Fe phase sat- 
urated in carbon. We assume that after the segregation 
process responsible for the generation of the liquid layer, 
the concentration in this layer will vary from 25%at. at 
the solid-liquid interface to 45%at. at the liquid-gas in- 
terface. We claim that the Benard-Marangoni convection 
governs the carbon transport in the liquid layer. The in- 
stability can be driven by temperature or concentration 
gradients; our calculations will demonstrate that concen- 
tration effects dominate. 

Since C-SWNTs grow in a bundle from one side of the 
NP, we neglect curvature effects, and work with a pla- 
nar geometry. We consider a liquid layer of thickness d 
and of infinite horizontal extent, with the same carbon 
concentration gradient as above. The free surface will be 
considered to be undeformable. Gravity effects can be 



neglected because Rayleigh-Benard convection is negli- 
gible compared to the Benard-Marangoni convection, as 
will be shown by our calculation. We model the Fe-C 
solution as a Newtonian and incompressible fluid. 

Since the BMI is driven by variations in surface tension 
a, which in turn are caused by gradients in concentration 
c, it is crucial to establish a mathematical relationship 
between these two quantities. Unfortunately, there is no 
experimental data available for such non-equilibrium sys- 
tems. Consequently, we extend the linear dependence 
fj = (To + 7c(c — Co) |l(Tj . valid in the eutectic region, to 
the entire hyper-eutectic region (c > 17.1%erf.), where 
carbon acts as a surface tension increasing solute. Here 
ao and Co are the eutectic composition surface tension 
and concentration, respectively. The solutal surface ten- 
sion coefficient, 7 C , is assumed constant. We assume the 
liquid layer is initially in an unperturbed state. The fluid 
is at rest and the initial concentration profile c is con- 
sidered linear with respect to the vertical distance Z, 
c = csl + 0Z, where (3 is a (positive) constant concen- 
tration gradient and csl is the carbon concentration at 
the solid-liquid interface. 

The motion of the liquid layer is described by the non- 
dimensionalized continuity and Navier-Stokes equations: 

V-m = 0, So' 1 {dtU + u-Vu) = -Vp + V 2 u, (1) 

where (x,y,z) are the three Cartesian coordinates, u — 
(u,v,w) is the velocity, p the pressure, Sc = p/ pD^ the 
Schmidt number, p, the viscosity, p the iron-carbon so- 
lution density and Dl the carbon diffusivity in the liq- 
uid phase. Here, we have scaled distance, time, velocity, 
concentration and pressure by d, d 2 /Dl, D^/d, (3d and 
/i£>i/ci 2 , respectively. We couple the above equations to 
a solutal diffusion process in the liquid layer modeled by: 

d t c + u-\7c = V 2 c-w, (2) 

where c = c — c is the deviation from the unperturbed- 
state concentration c. 

At the solid- liquid interface (z = 0), the velocity is 
zero, u = 0, there is no flux of carbon, d z c\ z =o = 0, 
and the carbon concentration is csl = 25%at. At the 
liquid-gas interface [z = 1), the surface is undeformable, 
u>| 2= i = 0, and the initial surface concentration of car- 
bon is clg = 45%at. Depending on the mass transport 
generated by the solutal BMI, c^g will change with time. 

At the onset of the instability, the boundary condition 
expressing the diffusion of carbon across the surface is 
given by d z c\ z= \ = 0. As the liquid layer is supersatu- 
rated with carbon at the onset of the instability, neither 
absorption nor adsorption are possible. Thus, the solutal 
Biot number Bi (characterizing rate of absorption) and 
the adsorption number Na [111] are both zero. 
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At the liquid-gas interface, the tangential stress is de- 
termined by the surface tension and takes the following 
nondimensional form: 

d zz w = -Ma s (d xx c + d yy c), (3) 

where the dimensionless solutal Marangoni number 

lcd{c L G ~ c L s) 



Ma. 
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(4) 



is an important predictor of the onset of the instability. 
The critical value of the solutal Marangoni number, Ma% 
above which the instability is generated in the case where 
Na and Bi are zero, is found from previous linear stability 
analyses 0,^3 to be Ma c s = 50. We will now show that 
Ma s > Ma 1 ^ under the particular synthesis conditions of 
C-SWNTs. 

We know from experimental observations that the di- 
ameter of the core of a NP is generally between 5nm 
and 20nm |4j . We postulate that the typical thickness of 
the liquid layer surrounding the core is between 2nm and 
4nm. The value of the solutal surface tension coefficient 
7 C has been determined in the literature 01 to be 0.03 
N/m ■ %at. We assume that the carbon diffusivity Dl 
in a highly supersaturated solution is equivalent to sub- 
stitutional diffusion because interstitial sites are mostly 
filled. In this case, the diffusivity is almost independent 
of the nature of the solute for liquid metals and varies 
between 10~ 9 m 2 /s and 10 _8 m 2 /s 0|. We choose the 
lower value of diffusion since the liquid layer is super- 
saturated. Finally, we use a viscosity within the range 
(i — 8 ± 2mPa ■ s, based on different extrapolations from 
data |l4j |. The solutal Marangoni number describing our 
system is then estimated to be between Ma s — 120 to 
400, which is above the critical Marangoni number Ma c s 
required to generate the instability. In contrast, the ther- 
mal and solutal Rayleigh numbers and the Marangoni 
thermocapillary number are estimated to be of the order 
of 10~ 20 , 10~ n and 10~ 3 respectively, far below their re- 
spective critical values. It follows from these calculations 
that the solutal Marangoni effect dominates all the other 
effects. 

Having described the onset of the Benard-Marangoni 
instability, we now determine the shape of the convec- 
tion cells. For this, we rely on the results of the weakly 
nonlinear analysis done by Bragard et al. The parameter 
necessary to determine the pattern is the distance from 
the threshold of the instability |10| : 



Ma s - Ma c s 
Mai 



(5) 



The value of e varies between 1.4 and 7 depending on the 
thickness of the liquid layer, and can only be considered 
to be approximate, given the imprecision in the experi- 
mental parameters. Using a linear extrapolation (Fig. 11 



and Table 1 from 0|), we determined the following sta- 
ble configurations for the convection cells when Na = Bi 
= 0: 



Stable configurations 


e 


Hexagons 


< e < 2.4 


Hexagons, rolls 


2.4 < e < 6.3 


Rolls, hybrid cells 


e > 6.3 



We conclude that the hexagonal stable configuration is 
favoured in the major part of the range of e. Note that 
we expect rolls and hybrid cells to be unfavorable for the 
nucleation of C-SWNTs. 

We now focus on the nucleation and growth of a C- 
SWNT at the center of a hexagonal convection cell. 
There is no mass transport between cells, and therefore 
it suffices to consider one convection cell. We know from 
fluid dynamics that in each convection cell, the liquid 
rises in the center of the hexagon and descends along the 
edges of the hexagon, forming a vortex ring 0] (Fig. 
The center of the top face and bottom face of a cell are 
stagnation points, where the velocity of the fluid is zero. 
The stagnation point at the bottom of a cell, lying on 
the solid-liquid interface, is of crucial importance since 
we postulate that a C-SWNT nucleates at this point. 
The tube grows by incorporating carbon atoms from the 
liquid. 

We envision the beginning of the nucleation process as 
follows: a fraction of the carbon atoms are transported 
close to the stagnation point on the carbide layer by the 
convective fluid flow (Fig|2Ji). The velocity of carbon 
atoms decreases close to the stagnation point, reducing 
their kinetic energy to below the adhesion energy on the 
metal catalyst. There is thus a semi-spherical influence 
zone around the stagnation point in which carbon atoms 
move slow enough to be adsorbed onto the catalyst. The 
semi-spherical shape is imposed by the axial-symmetry of 
the flow near the center of the cell. This influence zone 
acts as a heterogeneous nucleation site. Once the carbon 
is adsorbed at the solid-liquid interface within this zone, 
it crystallizes under its minimum energy configuration. 
It has been shown 16] that the nucleation of a closed 
hemispherical cap is favoured in presence of a metal sur- 
face. 

As the cap is formed, the topography of the bottom 
of the cell changes. The stagnation point is lost while 
a stagnation line (circumference) at the base of the cap 
is created (Fig. Eb)- The shape of the influence zone 
also changes: it becomes a sector of a torus surround- 
ing the stagnation line (Fig. |2>c). The hemispherical 
shape of the cap, observed in experiments 0], guarantees 
that the angle of contact between the cap and the carbide 
layer is nearly vertical. Consequently, the carbon atoms 
decelerate sufficiently and can be incorporated between 
the cap and the solid surface. The binding of the cap 
with the metal surface prevents its closure at the bottom 
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FIG. 2: Steps of the growth: a) Flow in one hexagonal convection cell, b) Formation of a C-SWNT cap on the solid-liquid 
interface, c) Growth of a C-SWNT in one hexagonal convection cell, d) Collective growth: a bundle of C-SWNTs 



and allows the root-growth [2| of a cylindrical structure, 
the C-SWNT (Fig. Et). We postulate that the carbon 
supply at the base of the cap, which is enhanced by the 
convective flow, exceeds the rate at which bonds in the 
cap can rearrange. This prevents the lateral growth of 
the cap and allows the C-SWNT extrusion with a con- 
stant diameter [J^. Since the size of a convection cell 
is known to be of the order of the liquid layer thick- 
ness [l5j, we can estimate the diameter of the C-SWNT 
to be a few nanometers. This agrees with the diameter 
of experimentally observed C-SWNTs, (0.7nm-3nm), 
Recent total energy calculations using density functional 
theory predict a larger diameter for the nanotubes.and 
conjecture a kinetic factor limiting this diameter, |l6| . 
Our BMI model provides precisely such a kinetic expla- 
nation. Note that outside the influence zone, it is more 
favourable for carbon atoms to follow the convection flow 
than to bind laterally on the C-SWNT since graphitiza- 
tion can only occur at much higher temperatures. 

The collection of hexagonal convection cells is respon- 
sible for the triangular lattice structure of the bundles of 
C-SWNTs (FIG. EJi) observed in experiments at the on- 
set of the growth Q. After the C-SWNTS grow beyond 
a certain point their organization may be influenced by 
other factors, but this stage of the growth is outside the 
scope of our model. Our model explains the nucleation 
and the beginning of the growth of C-SWNTs in the ab- 
sence of absorption (Bi = 0). However, the carbon atoms 
initially contained in the liquid layer are not sufficient 
for the growth of long bundles of C-SWNTs. Carbon 
atoms must therefore be absorbed from the surrounding 
gas to maintain the growth process. The extension of the 
present model to the case when Bi ^ is not trivial and 
will be addressed in a future study. 

In summary, using the Fe-C phase diagram, we have 
explained why a high cooling rate is necessary in the 
synthesis process to generate a layer of liquid metal sur- 



rounding the NP. Leveraging a linear stability analysis, 
we have shown that the solutal BMI can be generated 
in the liquid layer. Exploiting a weakly nonlinear anal- 
ysis, we have shown that the hexagonal pattern can be 
favoured under the synthesis conditions considered. Once 
initiated, the BMI is seen to be responsible for the nu- 
cleation and growth of C-SWNTs. 
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